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Abstract 

"*^ ■ A model is presented to estimate the fraction of Supernova Type-II events 

U ■ (SNII) occurring inside vs. outside a spiral arm for a given star formation 

^ ! episode. The probability distribution function (PDF) for this fraction is given 



»i^' for use in models similar to those of Shaviv et al. [13 11 . The calculated 



O I PDF for the SNII fraction, SNIIin/totai, defined as the number of SNII inside 

^ • a spiral arm divided by the total number of SNII from a star formation 

ci^ - event, provides a constraint on the magnitude of supernova remnant (SNR) 

concentrations used in cosmic ray propagation models attempting to explain 

K^ , the PAMELA anomaly. Despite the concentration of star formation within 

\Q ' spiral arms, this model predicts the majority of SNII events actually occur in 

On . inter-arm regions and calls into question the SNR concentration assumption 

of Shaviv et al. 

t^^ ■ Keywords: PAMELA Anomaly, Supernova Remnant 

o ■ 

1. Introduction 

" '"^ ' 

^ . Recently, the PAMELA collaboration released their results for the Galac- 

H i tic cosmic ray (GCR) positron fraction, J^+/{J^- + Je+)[1]. The remark- 

able feature in the PAMELA results is an increase in the positron fraction 
at energies above ~ 10 Gev. Standard GCR propagation models assume 
all positrons are secondarily produced via interactions with the interstellar 
medium (ISM). This results in a decreasing positron fraction with energy. 

Many models have been proposed to account for this discrepancy |^: DM 
particle annihilations, reacceleration of secondaries in supernova remnants 
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(SNRs) (eg. |2l]), inhomogeneities in the distribution of SNRs (eg. |13||ll[|). 
For a more thorough discussion of proposed models and their imphcations, see 
[GJ . Models relying on a concentration of SNR in the spiral arms of the Milky 
Way (MW) relative to the inter-arm regions have reproduced many observed 
results [ll[|l3[. However, the hterature is sparse on the specific value of this 
concentration; thus, it has been difficult to constrain the magnitude of the 
concentrations used to plausible ranges. 

In section 3, I present the model used to calculate SNIIin/totai as well 
as the necessary input parameters. Section 4 presents the results of the 
calculations and the PDF of SNIIin/totai- Section 5 provides discussion on 
the use of the results and possible refinements. 

2. Model and Inputs 

If we assume a dense cloud in the ISM collapses to form a cluster of stars 
upon its first contact with the density enhancement of a spiral arm we can 
compare the time needed for a newly formed star to orbit through and exit 
the arm to the lifetimes of stars in the cluster. The spread in stellar velocities 
and lifetimes naturally leads to some SNII events occurring within the spiral 
arm and some occurring after the progenitor has orbited out of the arm. 

Calculating the ratio of SNII events inside a spiral arm to the total number 
of SNII expected from a single star formation event, SNIIin/totai, is accom- 
plished by finding the ratio of two areas in the initial mass function (IMF) 
for star formation. The two areas are defined by 1., the mass above which 
stars explode in SNII, and 2., the mass at which stars explode via SNII while 
still within the spiral arm. This second mass-point on the IMF is a function 
of several parameters. 

Given that the pattern speed of the spiral density wave is lower than the 
orbital speeds of stars near the Solar Galactocentric radius, recently formed 
stars will, with time, exit the spiral arm along its leading edge (See Figure 
1). However, if a given star were massive enough it would evolve through 
its lifecycle and end in a SNII while still within the spiral arm. Thus, to 
calculate SNIIin/totai "we need to determine the fraction of stars with masses 
high enough to eventually explode as SNII that are also massive enough to 
do so prior to existing to confines of the spiral arm. 

Thus, calculating the SNII fraction requires knowledge of 3 observation- 



ally constrained parameters: 1) the shape of the IMF 10|; 2) the mass at 



which stars end their lives via SNII IJ]; 3) the mass at which stars go SNII 



within the spiral arm, which depends on Solar Galactocentric radius 12| . 



Spiral Pattern Speed J7(|,j5|], orbital speeds of newly formed stars [I4I, the 



size of the spiral arm |8|,|3|, and stellar lifetime as a function of stellar mass 



Each of these input parameters has an associated error. In order to ad- 
dress this, I evenly sampled the error range for each parameter and calculated 
SNIIin/totai for all possible combinations of sampled values. This maps out 
the changes in SNIIin/totai with respect to errors in all input parameters. The 
distribution of calculated SNIIin/totai values, when normalized to area=1.0, 
provides a probability distribution for SNIIin/totai- 

In addition to the ranges and iteration values given in Table 1, the con- 
version of stellar mass to stellar lifetime was performed for stellar metalicities 
of Z=[0.001], [0.02], [0.03] 0. 



Input Summary 


Parameter 


Range 


Iterations 


SNII Mass Limit, Mq 


7.0-9.5 


11 


High Mass IMF Slope 


2.3-2.7 


9 


Galoctcentric Solar Radius, kpc 


7.9-9.0 


12 


Solar Orbital Speed, km/s 


223.0-255.0 


9 


Spiral Pattern Speed, km/s/kpc 


5.0-18.0 


14 


Spiral Arm Size, kpc 


0.5-3.0 


14 



Table 1: Each input parameter was stepped evenly through its observational range to 
produce the listed number of iterations 



3. Results 

The PDF of SNIIin/totai is shown in Fig. 2. This distribution can be 
thought of as a measure of how far a given star travels from its birthplace 
in the spiral arm (in the arm's reference frame). Framing this distance in 
terms of the SNIIin/totai fraction folds the tails of the distribution into the 
final points at either end. 

A fit of the results, normalized for use as a PDF, is given in Eqn. 1. 
The distribution is fit with delta functions of finite area at the endpoints and 
between as an exponentiated polynomial. Table 2 contains values for the 
coefficients. 
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Figure 1: As a molecular cloud orbits into a spiral arm, gravitational perturbations induce 
star formation. Since the orbital speeds of the stars are greater than the spiral arm pattern 
speed, the stars can, depending on their lifespans, exit the arm's leading edge. The width 
of the spiral arm, L, and the difference in orbital speeds, AV, define a timescale for 
determining a given star's location at the end of its life. 
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aoe^'=i"'^* \Where / j(x) 



for x=0; 0.0903762 

for < X < 1; 0.8951162 
for x=l; 0.0145076 

(1) 



4. Discussion 

Models utilizing a concentration of SNII in the spiral arms to account for 
the PAMELA anomaly can be roughly tested by the integration of Eqn. 1 
for the range of concentration values which allow the model to fit observed 
data. This provides a simple constraint on the likelihood of the concentration 
required by the model. 

Piran et al. ll| assume a value of 4 times more SNII in spiral arms 
than out (4:1). This corresponds to a SNIIin/totai fraction as used in this 
paper of 0.8. The robustness of their model was not discussed with respect 
to this assumption, but if we assume a range of 3:1 to 5:1 still reproduces 
observations within their model then we can assess the likelihood. 
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Figure 2: SNILin/totai fraction vs. Probability. The fit used to produce the PDF is plotted 
over the points. Discontinuities occur at values of 0.0 and 1.0 (see Results section). The 
probability for 0.0 is 0.0904 



A concentration between 3:1 and 5:1 (| < SNin/totai 



< |) corresponds to 
0.029. This shows that a ratio of 4 times more 



an area under the PDF of 

SNII in a spiral arm than out is, in this model, not a robust assumption. 

Further, this illustrates a general prediction of this SNII model. Namely, 
for a significant fraction of the parameter space, there will, in fact, be a 
deficiency of SNII in spiral arms; most SNII will occur in inter-arm regions. 
The median calculated value is ~ 0.329, corresponding to roughly 2 SNII in 
the inter-arm region for every 1 SNII inside a spiral arm. 

The inputs used in this model are well constrained observationally, with 
the exception of spiral arm size. The physical extent of spiral arms is inher- 
ently variable; further, the edges are not sharply dehneated. Coupled with 
possible eccentricities in the orbits of stars and molecular clouds, the range 
of values for the effective size of the spiral arm used in the calculations rep- 
resents the only significant source of systematic error. Clearly, in the model 
smaller arm sizes directly bias results toward lower SNIIin/totai fractions with 
larger arm sizes tending in the opposite direction. 

In this paper, the large range in spiral arm size (see Table 1) was used 



Fit Parameters 


Parameter 


Value 


ao 




1.167706502 


ai 




-1.23164 


a2 




16.0944 


as 




-42.1807 


a4 




34.2674 


as 




-9.25471 



Table 2: 
to account for these issues. However, with hmits on the overdensity in the 



ISM required by the models of PiranpjJ and Shaviv[l3|, a better definition 



of spiral arm size as applied to calculating SNIIin/totai could be constructed. 
At that point, variances in orbital eccentricities at the Solar Galactocen- 
tric Radius, though relatively small, could be more readily accounted for to 
construct a more locally applicable version of this general model. 

With the majority of SNII occurring outside of the progentior star's par- 
ent spiral arm, models attempting to fit the PAMELA anomaly with spatial 
concentrations of SNII (rather than known locations of individual SNII rem- 
nants) must take note of the level of concentration needed and it's likelihood. 
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